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The Lagrangian statistics of heavy particles and of fluid tracers transported by a fully developed 
turbulent flow are investigated by means of high resolution direct numerical simulations. The 
Lagrangian velocity structure functions are measured in a time range spanning about three decades, 
from a tenth of the Kolmogorov timescale, r,,, up to a few large-scale eddy turnover time. Strong 
evidence is obtained that fluid tracers statistics are contaminated in the time range r G [f : fO]r,, by 
a bottleneck efl'ect due to vortex trapping. This effect is found to be signiflcantly reduced for heavy 
particles which are expelled from vortices by inertia. These findings help in clarifying the results 
of a recent study by H. Xu et al. [Phys. Rev. Lett. 96 024503 (2006)], where differences between 
experimental and numerical results on scaling properties of fluid tracers were reported. 
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Suspensions of dust, droplets, bubbles, and other finite- 
size particles advected by incompressible turbulent flows 
are commonly encountered in many natural phenomena 
and applied processes ranging from cloud formation to 
industrial mixing 0, Q. Understanding their statisti- 
cal properties is thus of primary importance. From a 
theoretical point of view, the problem is more compli- 
cated than in the case of fluid tracers, point-like parti- 
cles with the same density of the carrier fluid: inertia 
is responsible for the appearance of correlations between 
the particle positions and the structure of the underly- 
ing flow. It is well known indeed that heavy particles 
are expelled from vortical structures, while light parti- 
cles tend to concentrate in their core. This results in 
the appearance of strong inhomogeneities in the parti- 
cle spatial distribution, often dubbed preferential con- 
centration |3, y. This phenomenon has recentlygath- 
ered much attention both from a theoretical [5|, |6j and 
a numerical [J, |3| point of view. Progresses in the sta- 
tistical characterization of particle aggregates have been 
achieved by studying particles evolving in stochastic flows 
[5|, la, Isl and in two dimensional turbulent flows [lOj. 
Concerning single particle statistics, there has been con- 
siderable less attention for heavy particles [ll|, u^, in 
contrast with the numerous studies devoted to tracers 
[11 H H H m H E US 113 For example, it is 
known that Lagrangian velocity structure functions of 
tracers display high intermittent statistics [l^ [V^ [iSj, 
while there arc no results concerning heavy particles. 
The aim of this Letter is to compare the Lagrangian 
statistics of heavy particles with that of fluid tracers 



evolving in the same turbulent flow, by means of Direct 
Numerical Simulations (DNS). In particular, we shall fo- 
cus on the Lagrangian structure function of second or- 
der, important for stochastic modelling of particle tra- 
jectories [23, I23, and of higher orders, to account for 
intermittency effects. For Lagrangian studies, DNS offer 
the possibility to reach Reynolds numbers comparable to 
experiments, with a full control of both temporal and 
spatial properties of a very large number of tracers and 
heavy particles. 

We consider particles with a mass density pp much larger 
than the density pf of the carrier fluid. In this limit, par- 
ticles evolve according to |2J| 
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X{t) denotes the particle position, V{t) its velocity, 
u{x,t) is the fluid velocity and Tg = 2ppa^ /{Qpfiy) is 
the particle response time, where a is the radius and v 
the fluid viscosity. The Stokes number, which quantifies 
the degree of inertia, is defined as St ~ Ts/tjj, where 
Tri = {v/eY^^ is the eddy turnover time associated to the 
Kolmogorov scale and e the av erag e rate of energy injec- 
tion. Equation (^ is derived in |2J| under the assumption 
of very dilute suspensions, where particle-particle inter- 
actions (e.g. collisions) and hydrodynamic coupling can 
be neglected. Tracers correspond to the limit Tg -^ 0, i.e. 
to the evolution dX/dt = u{X{t), t). 
We performed a series of DNS of homogeneous and 
isotropic turbulence on a cubic grid with resolutions up 
to 512"^ (reaching a Taylor micro-scale Reynolds number 



R\ ~ 180) and transporting millions of tracers and par- 
ticles with St G [0.16:1]. Particles are initially seeded 
homogeneously in space with velocities equal to the lo- 
cal fluid velocity, already in a stationary configuration. 
Then they evolve according to equation ^ for about 2 
to 3 large-scale eddy turnover times before reaching a 
Lagrangian statistical steady state. It is only once the 
particle dynamics has completely relaxed that measure- 
ments are started. Details on the simulation parameters 
can be found in previous reports |llL l25l | . For fluid trac- 
ers we also present results obtained with resolution 1024"^ 
and R\ « 300; these data are described in |l7ll21|. 
The Lagrangian Structure Functions (LSF) 



S^P\T) = {[V{t + r)-V{t)r). 



(2) 



measure the time variations of any component V of the 
tracer or particle velocity. For time lags r in the incrtial 
range (i.e. when r,, ^ r <IC T^, where Tl is the integral 
Lagrangian time scale), the LSF of tracers are expected 
to display power law behaviors S^p\t) oc t^^p^ (see |lfil | 
and references therein). A dimensional estimate derived 
from the Kolmogorov 1941 theory for Eulerian turbulence 
predicts S^p'> (r) a (er)^/^ . Considerable deviations from 
the non-intermittent scaling C(p) = ip/2 were observed in 
both experiments [Ullil and DNS [llllll, with however 
some disagreement on the actual values of the exponents. 
For heavy particles there is not any reference theory. As 
the effect of inertia on the velocity statistics becomes less 
and less important when increasing t/ts |5J , we can guess 
that particles recover the statistical properties of tracers 
when T,, ~ Ts <C T <C Tl- 

We report in the sequel numerical measurements of 
the LSF of both particles and tracers for time lags 
ranging from r^/lO up to IOOt^ of the order of the 
integral time scale T^. Here we anticipate the main 
findings. Comparing the statistical properties of tracers 
and particles, we find further evidence that the formers 
are strongly affected by vortex trapping well above the 
Kolmogorov time scale, i.e. in the range r g [1:10]t^, 
while trapping becomes less and less important at in- 
creasing St. Trapping events spoil the scaling properties 
for time lags expected to be inside the inertial range of 
turbulence. This observation explains the disagreement, 
pointed out in |l8l | , existing between the scaling behavior 
of fluid tracers LSF measured in experiments [IJ, [ig 
and DNS [13, l21|. We argue that this discrepancy 
stems from the fact that the experimental LSF scaling 
exponents are measured in the time range [3 : 6]r^, 
exactly where trapping is effective. This leads to a more 
intermittent statistics than that measured in the inertial 

range [13 EH. 

Figure ^ summarizes the results for the second order 
LSF at varying Stokes. It should be noted that it is 
very difficult to identify a power-law scaling range for 
any Stokes and any time lags. This is evident from the 
absence of plateaux in the logarithmic derivatives of the 




FIG. 1: S^'^'{t) vs r/Trj in log-log scale for tracers (■) and 
heavy particles with St = 0.16, 0.37, 0.59, 1.01, and 1.34 
(solid lines from top to bottom). Bottom inset: compensated 
plot S"-'(r)/r vs r/rr,, same symbols as in the body. Top 
inset: logarithmic derivative dlog(5'^'(r))/dlog(T) vs r/rri, 
same symbols as in the body. For each Stokes number, aver- 
ages in ||2J are performed over N = 5 x 10^ trajectories that 
last for about 200t,,. All curves are obtained by averaging 
over the three components of the velocity vector to increase 
the statistics. 
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FIG. 2: Log-log plot of S^^^r) vs S^^\r) of order p = 4 and 
p — 6 for fluid tracers (■) and heavy particles with St = 0.16, 
0.37, 0.59, 1.01, and 1.34 (solid lines from top to bottom). The 
two straight lines represent the dimensional non-intermittent 
scaling (^(p)/^(2) = p/2. Notice the similarity of heavy par- 
ticles for different Stokes and the marked difference of the 
tracers with respect to particles for values close to the vis- 
cous scale, S^'^'{Tri) ~ 0.1 (vertical dotted line). 



LSF plotted in the top inset. The apparent plateaux in 
the bottom inset, where we show the compensated LSF 
normalized to the S^'^\t)/t at changing Stokes, should 
thus not be considered as an evidence of a linear power 
law. Still in the bottom inset, it is worth noticing 
the strong effects induced by inertia: already for the 
smallest Stokes number we considered {St = 0.16), we 
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FIG. 3: (Color online) Logarithmic derivative of the ESS plots of Fig. |21vs r/rn for p = 4 (left) and p = 6 (right). Data for 
fluid tracers at Rex ~ 180 and Re^ — 300 are plotted as ■ and •, respectively. The solid lines refer to heavy particles for 
St = 0.16, 0.37, 0.59, 1.01 (from bottom to top). The dotted straight line marks the Lagrangian Multifractal prediction obtained 
for tracers in [iTl (C(4)/C(2) — 1.71,(^(6)/("(2) = 2.26. The grey (yellow online) band corresponds to the experimental values 
measured for tracers in llal at comparable Reynolds numbers, and fitted in the range r £ [2 : 5]r^ (shown with the arrows). 
Errors have been estimated form the variations among the three velocity components. 



observe a clear departure from tracer behavior. 
It is now interesting to look for non-trivial effects at 
higher-order moments. As customary in turbulence 
studies, assessing deviations from a simple dimensional 
scaling can be done by comparing all moments against 
a reference one, a procedure originally proposed for 
Eulerian structure functions and dubbed Extended Self 
Similarity (ESS) [26|. This amounts to study the scaling 
behavior of the order p LSF as a function of that of 
order 2. used as a reference for Lagrangian statistics. 
This procedure is known to decrease finite-Reynolds 
effects and is frequently used for measuring scaling 
exponents in experiments and simulations. The price 
to pay is that only relative scaling exponents C(p)/C(2) 
can be measured. Figure |21 shows S'^p^t) as a function 
of S'^-^^t) for p = 4 and p = 6 and various values of 
St. Two observations can be done. First, the LSF of 
both inertial particles and tracers have an inertial-range 
scaling behavior that deviates significantly from the 
dimensional one C,{p)/C,{2) — p/2. Second, the differ- 
ences between tracers and particles around t ^ t^j, 
i.e. for S^'^^Tr)) ~ 0.1, are now even more pronounced. 
To be more quantitative in Fig. the local slopes 
dlog(S'(P)(r))/dlog(S'(2)(T)) of all curves are shown. Let 
us first focus on the behavior well inside the inertial 
range. For IOt,, < r < T^, all the data sets display a 
tendency to converge, within error bars, to the same 
scaling behavior as that of fluid tracers, irrespectively of 
the value of the Stokes number. The large error bars are 
due to unavoidable large-scale anisotropic fluctuations 
in the statistics. In the same figure we also report 
the multifractal prediction for fluid tracers obtained 
in [13 (see also |l3ll27|). by translating the well-known 
Eulerian multifractal model to the Lagrangian domain. 



Notice that the inertial-range behavior is compatible 
with the multifractal prediction also for heavy particles. 
For time lags in the range t^i < t < lOr^, the local 
scaling exponents reveal the presence a strong bottleneck 
effect, much more pronounced for tracers than particles. 
At these time scales, both for p = 4 and p = 6, fluid 
tracers have very intermittent scaling properties that 
are however smoothed out as soon as inertia is switched 
on. This is due to the fact that tracers {St = 0) 
may experience vortex trapping lasting for rather long 
times, while inertial particles (even for S'i <C 1) are 
expelled from vortex filaments. The values for the 
scaling exponents of the LSF given in the experimental 
study of fluid tracers in [18j at comparable Reynolds 
numbers are also represented in Fig. |31 in the shaded 
band. In Ref. [l8J, the scaling exponents are measured 
in the range r G [2 : 5]r^ (see arrows in the figure), 
where trapping into vortex filaments gives the dominant 
contribution leading to a more intermittent statistics. 
For these time lags, the values of the scaling exponents 
reported in [ig are in good agreement with DNS results. 
However, these exponents substantially differ from the 
inertial-range values of the logarithmic derivative that 
we observe at r > lOr^. 

It is worth stressing that the statistical signature of par- 
ticle trapping into vortex filaments has already been the 
subject of experimental and numerical studies of fully 
developed turbulent flows. In particular, it was shown 
that such events play a crucial role in determi ning the 
intense fluctuations of tracer acceleration [l4L Il6l llTl l28l| . 
This effect has previously been highlighted by filtering 
out the intense vortical events in the LSF (see Fig. 3 of 
[21| and compare it with Figs.na-ndO- Such a filtering 
is here obtained dynamically by switching on inertia, 



whose major landmark at small values of the Stokes 
number is to expel heavy particles from vortex filaments. 
In conclusion, we have shown that for large-enough time 
lags (r > lOr^) the scaling properties of heavy particles 
with St < 1 tend to be almost independent of St, 
within error bars. The relative scaling exponents of LSF 
for the investigated Stokes numbers are found in the 
range C(4)/C(2) £ [1.55 : 1.75] for the fourth order and 
C(6)/C(2) G [1-8 : 2.4] for the sixth one. Investigations 
at larger Re\ are needed to confirm the robustness of 
the universality enjoyed by heavy particles for large 
time lags. This observation can be explained by the fact 
that the effect of inertia on the particle velocity should 
disappear progressively when increasing the scale [3 ■ For 
time lags well inside the inertial range and much larger 
than the response time Ts , the delay of the particles on 
the fluid motion becomes negligible and a "tracer-like" 
physics is recovered [23. For the case of fluid tracers, 
it is important to be very careful in assessing scaling 
properties because trapping into vortex filaments may 
spoil the inertial range scaling behavior for time lags up 
to r ~ lOr^ and even more. Such trapping effects arc 
less effective for heavy particles due to their dynamical 
properties that bring them outside of strong vortex 
filaments. These results are important for stochastic 
modelling of both tracers and heavy particles, where the 
statistical properties of velocity and acceleration along 
the trajectories are the main ingredients for developing 
appropriate models. 
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